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Abstract 

General point interactions for the second derivative operator in one di- 
mension are studied. In particular, PT-self-adjoint point interactions 
with the support at the origin and at points ±Z are considered. The 
spectrum of such non-Hermitian operators is investigated and condi- 
tions when the spectrum is pure real are presented. The results are 
compared with those for standard self-adjoint point interactions. 
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1 Introduction. 



Exactly solvable models are used in quantum mechanics to obtain Hamiltonians describ- 
ing realistic physical systems but having the important property of being exactly solv- 
able, i.e. that all eigenfunctions, spectrum and scattering matrix can be calculated in 
closed form using elementary functions. A large class of such operators can be obtained 
using the method of point interactions described in detail in 0] in application to the 
theory of self-adjoint operators. It is not surprising that point interactions can be used 
to obtain exactly solvable non self-adjoint operators, see for example where the 

Schrodinger operator on the half-axis with complex boundary condition at the origin has 
been considered as a model for dissipative operators. Using dissipative operators instead 
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of self-adjoint ones one obtains Hamiltonians describing irreversible quantum systems 
which are of great importance especially for evolution problems. It has been discovered 
recently that Hamiltonians possessing so-called PT-symmetry can be non self-adjoint and 
have the real spectrum at the same time. This discovery observed first using numerical 
computations has been approved by extensive analytical calculations. Many Schrodinger 
operators which have been studied have complex PT-symmetric potentials but real dis- 
crete spectrum. Even if no proof that the spectrum is real for the whole class of operators 
has been discovered, it was expected that the relations between these two properties are 
rather close. In the present article we consider T^T-symmetric operators with point inter- 
actions, which are exactly solvable. This gives us an opportunity to study the relations 
between the PT-symmetry of the Hamiltonian and the reality of its spectrum in full de- 
tail. We show in particular that none of these two properties imply the other one, but 
that there exist Hamiltonians with point interactions which are PT-symmetric and have 
real spectrum. Since the operators under investigation are not self-adjoint the classical 
von Neumann theory cannot be applied without modification. Instead the method of 
boundary conditions is used. It is shown how the notion of PT-symmetry has to be mod- 
ified in application to the extension theory for linear operators. In addition we introduce 
the notion of PT-self- adjoint operators. The operators obtained can be used to describe 
irreversible systems in the framework of the recently appeared PT-symmetric quantum 
mechanics described below. 

Let us start by presenting recent results concerning one-dimensional Schrodinger op- 
erator with PT-symmetric potentials determined by the following expression 

L=-^ + KW, (1) 

where the potential V{x) is not a real valued function, but satisfies: 



V{-x) = V{x). 

Using the space parity operator V : {Vip){x) — ip{—x) and the complex conjugation 
operator Tip — the last property can be written as 

-prv = V. 

Then the differential operator L is formally PT-symmetric 

VT L^L VT. 

Interest to differential operators with such symmetry became enormous, after it has been 
discovered that some of these operators have real spectrum (like self-adjoint operators) 
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and therefore can be used in construction of a new, PT-symmetric quantum mechanics.^ 
This fact was first discovered numerically by CM. Bender, S.Boettcher and collaborators 
11; H; S B 0' mi]- Analytical studies of PT-symmetric Hamiltonians were carried out 
by M.Znojil and his collaborators |2l|, |2^, Numerous non-self- 



adjoint operators with real spectrum were studied (see e.g. [|^, [H, [T^, |TI 

The aim of the current paper is to describe exactly solvable "PT-Hermitian operators 
constructing using the method of point interactions 0, Standard symmetries of (self- 
adjoint) point interactions in dimension 1 were studied in We are going now to 
extend these results to include PT-self-adjoint point interactions. Starting from the 
second derivative operator in L2{M) (which is PT-self-adjoint) we construct non trivial 
point interactions leading to operators having "PT-symmetry. Their spectral properties 
are described in Section 3. Section 4 is devoted to construction of local point interactions 
concentrated at different points on the real line. One of the main questions investigated 
is the relations between the PT-symmetry of the operators and reality of their spectrum. 
The family of point interactions with real spectrum is characterized in Section 5. 



2 PT-self-adjoint point interactions. 



The method of point interactions is well described in several monographs This 



method is based on the fact that differential expression which is formally symmetric (or 
PT-symmetric in our case) does not determine the operator in the Hilbert space uniquely. 
To determine the operator one has to specify its domain. One can start from a certain 
standard differential operator (i.e. Laplace operator) and restrict it to a certain densely 
defined operator. Then extending the restricted operator to another PT-Hermitian op- 
erator one can get an operator with nontrivial spectral structure. Thus the main tool to 
be used is the extension theory for linear operators. The method we are going to present 
is similar to the method of point interaction developed for self-adjoint operators. There- 
fore let us introduce two definitions which are similar to the definition of symmetric and 
self-adjoint operators in conventional quantum mechanics. 



Definition 1 An operator L is called PT-symmetric0 if it is densely defined and the 
following inclusion holds 

L* D VLV. (2) 

^It is not true in general that the spectrum of every T'T-symmetric operator is real. 

^The operators described by these definitions should be better called T'-symmetric and T'-self-adjoint, 
but we prefer to use notation T^T-symmetric and PT-self-adjoint in order to underline the entire relations 
with ■PT-symmetric quantum mechanics already described in the literature. 
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A VT -symmetric operator L is called PT-self-adjoint if and only if 

L* = VLV. (3) 

The second derivative operator L = — ^ with the standard domain W2{Si) is both 
self-adjoint and PT-self-adjoint. In the current section we study point perturbations at 
the origin of this operator leading to PT-self-adjoint operators. By point interaction we 
mean any interaction which vanishes on the functions with the support separated from 
the origin. In other words a linear operator A is a point perturbation at the origin of the 
operator L if and only if its restriction to C^{1R \ {0}) coincides with the restriction of 
the operator L: 

A.\c^{iR\{o}) = L|c^(jR\{o}) = Lq. 

The operator Lq determined by the last equation is both symmetric and PT-symmetric. 
Moreover 

VLoV = Lo, VLIV = LI 

and thus any PT-symmetric extension of Lq is a restriction of the maximal operator 
Lma.x = LI being the second derivative operator in L2{M) with the domain Dom(Linax) = 
W2{M \ {0}). All such PT-self-adjoint extensions and therefore all PT-self-adjoint point 
interactions at the origin are described by the following 



Theorem 1 The family of VT -self- adjoint second derivative operators with point inter- 
actions at the origin coincides with the set of restrictions of the second derivative operator 
Lma.x = defined originally on W2{JR \ {0}), to the domain of functions satisfying 

the boundary conditions at the origin of one of the following two types 

^ ( ^(-0) ^ ■ m 

I i^'i+o) i [ ^'(-0) i ' 



with the matrix B equal to 

B 



je 



VYTbce^'t' b 

c VI + be e-'^ 



with the real parameters b > 0,c> —1/6,0 0,(j) & [0, 2tt); 
II. 

/ W(+0) = /iie*V(+0) 



I ho^'i+0) = -hie-'y{-0) 

with the real phase parameter 9 G [0, 27r) and with the parameter h = (/iq, hi) taken from 
the (real) projective space P^. 

^If the parameter h is equal to zero, then the last inequality can be neglected. 
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Proof We have already proven that any PT-self- adjoint extension of Lq is a restriction 
of Lmax- Consider the vector space of boundary values of functions from the domain 
of Ljnax and the map F adjusting to any function ip e \ {0}) its boundary values 



r : -0 



/ V'(+o) \ 

^'(+0) 
^(-0) 



The fundamental symmetry operator V acting in the space of boundary values coincides 
with the operator of multiplication by the following matrix P 



PT = TV, where P 



( 

1 
V 









-1 



1 








\ 

-1 



) 



(6) 



Similarly the operator of complex conjugation T is mapped by F into the operator T of 
complex conjugation in £7^. 

The closure of the operator Lq is defined on the functions having trivial boundary 
values at the origin and Vice Versa any function from the domain of Lmax with triv- 
ial boundary values belongs to Dom(Lo). Thus the dimension of the quotient space 
Dom (Lmax)/Dom (Lq) is equal to 4. Any PT-self-adjoint extension of Lq can be de- 
scribed as the restriction of Lmax to the set of functions with boundary values from a 
certain 2-dimensional subspace £ of C^. Every such subspace can be described by two 
(linearly independent) boundary conditions using a certain 2x4 dimensional rank 2 matrix 
Q = {qij} as follows 



g2iV'(+0) + g22V''(+0) + g23V'(-0) + g24V''(-0) = 



(7) 



The restriction of the maximal operator to a two-dimensional subspace C possesses the 
T'T-symmetry if and only if 

X eC^ PTX e C. 



Let us prove now that the PT-self- adjoint point interactions can be parameterized in 
at least one of the following four ways 
A 

b>0,c> -l/b,e,(pe [0,27r); 



(8) 



^/T+bce''f' 



VI + be e 



■iff) 



i'i-0) 



(9) 
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6 > 0,c > 
C 

a,b> 0,6, 
D 



-1/6,^,0 G [0,27r); 



V^(-O) 



ie 



e [0,27r); 



ae 
-be 



ae 
-be~ 



be''^ 
— ae~' 



^'(+0) 



^(+0) 
^(-0) 



(10) 



fill 



a,b>0,9,^E [0,27r). 



Since the matrix Q appearing in (|^) has rank 2, at least one of its six 2 x 2 minors is 

non-degenerate. Depending on which minor is non-degenerate the boundary conditions 

a P \ 



(^ can be written in different ways using certain 2x2 matrices B 
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1) det 

\ q2i q22 



2) det 

3) det 

4) det 

5) det 

6) det 



Qn qi3 

Q21 Q23 

qn qi4 

q2i q2A 

qi2 qi3 

q22 q23 

qi2 qu 

q22 q2i 

qi3 qu 

q23 q2i 



7^0 

7^0 
7^0 



^-(+0) 
^'(+0) 
^(+0) 
^(-0) 
^(+0) 

^(-0) 
^'(+0) 

^'(-0) 
^(-0) 

V^'(-o) 



: B 

: B 

-- B 

-- B 

-- B 



■ ^(-0) 

V'(-o) 

^(+0) 

' ^(+0) 
^(-0) 

■ ^(+0) 



Our aim now is to characterize all matrices B leading to PT-symmetric boundary condi- 
tions. 

Consider the first case. Suppose that the function ip satisfies the boundary conditions. 
The boundary conditions for the function VTip are given by 



-a —13 
-7 —5 



^(+0) 
^'(+0) 



These conditions coincide with the original one if and only if 

-1 



—a —13 
—7 5 



a (3 
7 5 



(12) 



Investigating cases 3), 4) and 6) one arrives to exactly the same equation on the matrix 
B. 

Equation ( P^ implies that det B = aS — ^ and the following four equalities hold 
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detB"' 



6 



detS^^' 
1 



detS 



a. 
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These equations imply in particular that the determinant has absolute value 1, i.e. there 
exists 9 G [0, 2tt) such that det B = e^*^. Consider then the matrix 

with the unit determinant a'6' — j3''~f' = 1. The entries of B' satisfy the following equations 

a' = 6', = f3', 

f = y, 6' = a'. 

Hence all such matrices B' can be parameterized by 3 real parameters: 

positive number b; 

real number c, c > —1/6; 

phase parameter G [0, 27r); 
using the following formula 

This implies that the boundary conditions can be written in the form in the cases 1) 
and 6) and in the form (P) in the cases 3) and 4).^ 

It remains to study the cases 2) and 5). The boundary conditions determine a VT- 
symmetric operator if and only if the following equality holds in both cases 



7 5 i [Pa 



(14) 



This equality is satisfied if and only if 

j a = —S 
\ /3 = -7 ' 

and such matrices can be parameterized by the following parameters: 

two positive numbers a, 6 G iR+; 

two phases 6,(f) G [0, 27r); 
using the formula 



B 



ae 



We thus get the parameterizations (p!0| ) and (|TT]). 

It remains to prove that each of the boundary conditions B, C and D can be written 

as (^) or (|^). (The boundary conditions of type A coincide with (^).) The boundary 

''One has to take into account that det B ^ and therefore the boundary conditions 1) and 6) and 3) 
and 4) are pairwise equivalent. 
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conditions of type B can be written as (^) if 1 + 6c 7^ simply by expressing ip'{—0) as a 
function of ^'(+0) ^iid iIj'{+0) from the second boundary condition (|^). If 1 + 6c = then 
the boundary condition is of type (^). The cases C and D can be studied similarly. 

Thus we have proven that these boundary conditions describe all maximal PT-symmetric 
extensions of Lq. One can easily check that all these operators are PT-self-adjoint as well. 
The theorem is proven. Box 

Similar results in the theory of self-adjoint point interactions are usually proven using 
von Neumann extension theory [1^, 26 1. Instead of developing its counterpart for VT- 
symmetric operators we preferred to give a constructive proof using boundary conditions. 

In what follows the boundary conditions given by will be called connected, since 
these conditions connect the boundary values on the left and right hand sides of the origin 
of functions from the domain of the operator. The boundary conditions of the second 
type given by (H) will be called separated. These definitions and the main proposition 
are quite similar to corresponding description of self-adjoint point interactions given in 
T9(| . The operators appearing in the decomposition corresponding to separated boundary 



conditions are just the second derivative operators on the half lines with complex boundary 
condition of the third type at the origin. Such non-self-adjoint operators can easily be 



studied (see e.g. [p5|] ). 

Thus the set of PT-self-adjoint point interactions can be parameterized by 4 real pa- 
rameters. The phase parameter 6 is redundant in the sense that the operators correspond- 
ing to different values of this parameter are unitary equivalent .0 The parameterization 
used in (|D and (^ is not optimal in the sense that the correspondence between the pa- 
rameters and the boundary conditions is not one-to-one (for example, if 6 = c = 0, then 
changing the phases by vr we get the same boundary conditions, similar problem occurs 
for separated boundary conditions), but we prefer not to dwell on this point. Connected 
boundary conditions given by the matrix 

B = e'^(^^^ a,b,ce IR,e e[0,n),a^ -bc=l (15) 

determine operators which are both self-adjoint and PT-self-adjoint. Separated boundary 
conditions leading to such operators are given by formula (^ with 9 = 0. 



3 Spectral problems. 

In this section we are going to study the spectrum of second derivative operators with 
PT-self-adjoint point interactions at the origin. Our main result can be formulated as 
^Similar fact for self-adjoint point interactions is described in full details M, [T^. 
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follows 



Theorem 2 The spectrum of any VT -self-adjoint second derivative operator with point 
interactions at the origin consists of the branch [0, oo) of the absolutely continuous spec- 
trum and at most two (counting multiplicity) eigenvalues, which are real negative or are 
(complex) conjugated to each other. 



Proof. Let us denote by A any PT-self-adjoint second derivative operator with point 
interaction at the origin (described by Theorem 1). To prove the theorem we are going 
to calculate its resolvent. Consider the resolvent equation 



{A-\)U 



F e L2{]R). 



The unique solution to this equation is the function U from the domain Dom {A) of the 
operator A satisfying the differential equation 



U" -XU = F, 



(16) 



everywhere outside the origin. We denote by k the square root of the energy parameter 
A = fc^, determined uniquely by > 0. Let us introduce two functions 

Jkx 2; > , , f 0, X > 

e~ 



e+{x) 



e_(a;) 



0, X < ' ^"'"^^ 1 e-*'^^, X <0 



Then any solution (from L2{1R)) to (|lg) can be written in the form 

foo Qik\x-y\ 



U{x) 



-F{y)dy + p+(F)e+(a;) + p_(F)e_(x) 



(17) 



-^oo 2ik 

where p± are two parameters to be calculated. The boundary values of the function U 
are 



f/(+0) 
f/'(+0) 
U{-0) 
U'i-0) 



1 

2ik 



{f- + f+)+P+ 



2ik 



if- + /+) + P- 



-U:-f^)-tkp. 



where f± = Jj^^^ F{y)dy. 

Consider the case of separated boundary conditions. Then 

/ 1 -e*^ fVl + hce}^ - ikb] 

V 



implies 



ik — e*^ (^c — %k\/\ + bee 
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_i_ (_i + e^o^yYTVce^'^) {f- + /+) + e*'|(/- - /+ 

I ^''i-kU- + U) + \ (-1 + e^'y/TTVce-^^) (/_ - /, 
This system of equations on p± is solvable if and only if the determinant of the matrix at 
the left hand side is different from zero 

- e'^ (c - 2ik cos 0Vl + hc- k%) ^ 0. (18) 

Thus the resolvent equation can be solved for all nonreal k, which are not solutions to 
the following quadratic equation 

bk^ + 2i cos (j)Vl + bck-c = 0. (19) 

The two solutions to the last equation 

.COS0V1 + be , v/(l - cos2 (f))bc - cosV 

fcl o = -t ; ± 



b b 

are either pure imaginary or symmetric to each other with respect to the imaginary axis. 
The corresponding energy values are real, or conjugated to each other. Hence the domain 
C \ (iR+ U {A^i, A;!}) belongs to the regularity domain of the operator L. 

Formula ([T7|) shows that the difference between the resolvents of the operators A and 
the unperturbed second derivative operator L has rank two. The spectrum of the (self- 
adjoint) operator L is pure absolutely continuous and fills in the interval [0,oo). The 
complement to the absolutely continuous spectrum is simply connected. Therefore the 
perturbed operator A has the branch of absolutely continuous spectrum [0, oo) as well 
||2^ , Let us study the singularities of the resolvent corresponding to the numbers fc^2- 
Let ^ki > 0, then the function 



e'^ i^Vl + bce'^ - ikibj e*^i^, x > 0, 

is a square integrable solution to the equation —ip" = kfipi, x 7^ and satisfies the 
boundary conditions (§). Therefore this function is a (discrete spectrum) eigenfunction for 
the operator A. Similarly the function ip2 = VTipi is an eigenfunction corresponding to the 
eigenvalue ki . If ki is pure imaginary {kf negative real) then the function ipi can be chosen 
PT-symmetric or -antisymmetric. In the special case (1 — cos^ 0)6c — cos^ = 0, cos < 
the two solutions to the dispersion equation (|19D coincide. Nevertheless the corresponding 
eigenvalue has multiplicity 1, since the boundary conditions (|^) are connected. Therefore 
we conclude that every solution to the dispersion equation ( pj]) from the upper half plane 
Qk > determines a simple eigenvalue fc^ of the operator A. Solutions to (|I9]) lying on the 
nonphysical sheet Qk < do not determine any eigenvalue of L, since the corresponding 
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solutions to the differential equation do not belong to the Hilbert space. The theorem is 
proven for connected boundary conditions. 

The proof for separated conditions is quite similar. The only difference is that in this 
case the operator A can be presented as an orthogonal sum of two operators in L2{M±). 
Each of these operators can have one (complex) eigenvalue and these eigenvalues are 
conjugated to each other. Therefore if the operator A has a negative eigenvalue then this 
eigenvalue has always multiplicity 2. □ 

In particular we have proven the following 

Corollary 1 The eigenvalues corresponding to VT -symmetric eigenfunctions of the op- 
erator A are real and negative. Every eigenfunction corresponding to any real eigenvalue 
of the operator A can be chosen VT -symmetric or -antisymmetric. 

This proposition holds in fact for any PT-self-adjoint operator with non-singular in- 
teractions. Note that the spectrum of the operator A is not always pure real. Let us 
study the case when the spectrum is real in more detail. The spectrum is pure real if 
and only if the two solutions to equation (|19D are pure imaginary or are situated on the 
nonphysical sheet '^k < 0. The two solutions are pure imaginary only if the discriminant 
is negative bcsm"^ (f) < cos"^ (p. If the discriminant is positive {bcsiia^ (p > cos^0), then the 
imaginary part of the solutions is given by _£2i^^2±i£ ^nd is negative only if cos0 > 0. 
We have proven the following 

Proposition 1 The spectrum of the VT -self- adjoint second derivative operator with con- 
nected point interaction at the origin is pure real if and only if the parameters appearing 
in (0j satisfy in addition at least one of the following conditions 
I. 

be sin^ < cos^ 0; 

//. 

6csin^ > cos^ and cos0 > 0. 

The spectrum does not depend on the phase parameter 9. The operators corresponding 
to different values of 6 are unitary equivalent (see [0], where similar results are proven for 
self-adjoint point interactions). 
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4 Local PT-self-adjoint point interactions. 



We consider now PT-self-adjoint operators with point interactions at positions x = ±Z, 
I G iR+. Every such operator coincides with a certain restriction of the second derivative 
operator — defined originally in L2{1R) on the domain W2{JR \ to the set of 

functions satisfying some boundary conditions at the points x = ±/. We restrict our 
consideration to the case of local connected point interaction without aiming to describe 
all possible restrictions leading to PT-self-adjoint operators. Suppose that the functions 
from the domain of the restricted operator satisfy the following conditions at x = / 

(20) 



where the matrix B 






G GL(2, C). Then to make the operator PT-symmetric 



we have to suppose that the following boundary conditions are introduced aX x = —I 







-1 



B 







-1 



M/(-/+) 
^'(-/+) 



(21) 



The restriction of the second derivative operator to the set of functions satisfying condi- 
tions (pOf ) and ( pTD will be denoted by A in this section. Let us study the spectrum of 
this operator. The resolvent of the operator A can be calculated explicitly using methods 
of the previous section. One concludes that the difference between the resolvents of the 
operators A and L has rank four and therefore the operator A has a branch of absolutely 
continuous spectrum [0, oo). Let us study the discrete spectrum of A. To calculate the 
eigenfunction one can use the following Ansatz: 

-ik{x+l) rf. ^ ^. 



cie 



iIj{x) 



C2 COS k{x + /) + C3 sin k{x + I), —l<x<l] 



ik{x—l) 



(22) 



X > t 



where k"^ = A, ^k > and ci, C2, C3, C4 are arbitrary complex parameters to be determined. 

!l|) gives the following 



Substitution of this function into the boundary conditions (|20|) and 
linear system on Cj 



/ 


1 


—a 


(3k 






( ^ 




ik 


-7 


5k 














a cos 2kl — /3k sin 2kl 


a sin 2kl + /3k cos 2kl 


-1 




C3 


V 





7 cos 2kl — 6k sin 2kl 


7 sin 2kl + 6k cos 2kl 


-ik ) 




^ C4 ; 



0. 



This system has a nontrivial solution if and only if the determinant of the 4x4 matrix is 
equal to zero, which gives us the dispersion relation 

sm2kl {-k^'lPl^ - ik'^ipS + P6) + kWa\^ - \6\^) + ik{a^ + aj) - I7P} 



+k cos 2kl {k'^ {a/3 + a/3) + ik{a6 + a6 + /3'y + /3-f) - (75 + 75)} 



(23) 



0. 



12 



Each solution to this equation from the upper half plane determines an eigenvalue of the 
operator A. 

In what follows we are going to study the operator with the interaction given by a 
sum of two delta-functions given by the following formal expression 

- + (m + w)6{x -l) + {u- iv)6{x + I), (24) 

where u,v & M. This potential is formally PT-symmetric. This operator is determined 
by the boundary conditions chosen as follows^] 

a = 7 = 1, P = 0, 6 = u + iv. 

In the case t> = we get a self-adjoint operator. Let us study the case u = 0, when the 
potential is pure imaginary. The dispersion relation takes the following simple form 

sin2A;/ [k^{l - v^) + 2ik - l} = 0. 

All solutions to this equation from the open upper half plane are determined by the zeroes 
of the second order polynomial in the brackets 

ki,2 = -i-^— , V ^ ±1. 

These solutions are always pure imaginary. Depending on the absolute value of v one or 
none solutions are situated on the physical sheet > 0. It is interesting to observe that 
the spectrum of this operator is always pure real. 



5 Point interactions with real spectrum. 

In this section we are going to discuss connected boundary conditions at the origin leading 
to second derivative operators with real spectrum. Namely the operator A = — ^ with 
the domain 

Don,(.,^{,.M^,^\{0},;(|(^)^(; f)(j<:^>)} (25, 



will be studied. We suppose that the matrix B = 




appearing in the boundary 



conditions is non degenerate (from GL{2, C)). Again it is easy to prove that the operator 
has branch of absolutely continuous spectrum [0, oo) and it remains to study its discrete 
spectrum only. Let us use the following Ansatz for the eigenfunction 

^'-) - { : > " • > »■ f^^) 

®To define the operator correctly one can use for example its quadratic form. 
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corresponding to the energy A = fc^. Substituting this function into the boundary condi- 
tions we get the 2x2 hnear system 



a — ik/3 — 1 \ / ci 
7 — ikS —ik 1 \ C2 



0, 



and the dispersion equation 

k'^(3 + ik{a + 5) - 7 = 0. (27) 
The last equation has the following two solutions 



_ -i{a + 6)± ^-{a + 5)2 + 47/? 
^^'^ - 2P ' 

if /? 7^ 0. Solutions of the equation situated in the closed lower half plane Qk < do not 
determine any eigenfunction of the operator A. Hence the spectrum of the operator A is 
real if and only if both solutions /ci 2 to equation ( pTf ) satisfy at least one of the following 
two conditions: 

1. SA;i,2 < 0; 

2. sRfci,2 = 0. 

The set of coefficients a, jS, 7, 5 satisfying the first condition can be parameterized by 8 
real parameters and leads to operators A with pure absolutely continuous spectrum [0, 00). 
Pure imaginary solutions to ( P?] ) can lead to nontrivial discrete spectrum. Therefore let 
us study the set of coefficients satisfying the second condition in more details. 

The solutions ki^2 are pure imaginary if and only if the following conditions are satis- 
fied: 

J-^^- ^^^^ 

The first two conditions imply that the complex numbers r = a + 6, f3, and 7 have the 
same phase. Therefore let us introduce the following parameterization: 

r = te*^, p = be'\ 7 = ce*^ 

where t, b, c are real numbers. In order to guarantee that solutions ki^2 are pure imaginary 
the real parameters should satisfy the inequality 

We conclude that the second family can be parameterized by 6 real parameters. The 
intersection between the two families of parameters leading to operators with pure real 
spectrum is not empty, but the second family is not included in the first one. 
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The four-parameter family of self-adjoint (connected) boundary conditions is contained 
in the second (6-parameter) family just described. Parameters leading to PT-self-adjoint 
boundary satisfy the first two conditions (|28|) . Hence the family of PT-self-adjoint (con- 
nected) boundary conditions leading to operators with real spectrum belongs to the second 
family of boundary conditions as well. 

6 Conclusions. 

The relations between PT-self-adjoint, self-adjoint and real second derivative operators 
with point interactions have been studied. One has to investigate whether the operators 
with real spectrum are similar to certain self-adjoint operators or not. It is easy to 
generalize these results in order to include point interactions with the support at arbitrary 
finite (or even infinite) number of points. The results obtained can be applied to construct 
exactly solvable few-body systems with unusual symmetries (following ^). This 
method can be generalized to include point interactions in spaces of higher dimension and 
higher order differential operators (like in [|l2l). 
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